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Measurements of forward proton production with incident protons and 
charged pions on nuclear targets at the CERN Proton Synchroton. 



Measurements of the double-differential proton production cross-section d^a / dpdQ in 
the range of momentum 0.5 GeV/c < p < 8.0 GeV/c and angle 0.05 rad < 6 < 0.25 rad 
in collisions of charged pions and protons on beryllium, carbon, aluminium, copper, tin, 
tantalum and lead are presented. The data were taken with the large acceptance HARP 
detector in the T9 beam line of the CERN Proton Synchrotron. Incident particles were 
identified by an elaborate system of beam detectors and impinged on a target of 5% of a 
nuclear interaction length. The tracking and identification of the produced particles was 

performed using the forward spectrometer of the HARP experiment. Results are 
obtained for the double-differential cross-sections mainly at four incident beam momenta 
(3 GeV/c, 5 GeV/c, 8 GeV/c and 12 GeV/c). Measurements are compared with 
predictions of the GEANT4 and MARS Monte Carlo generators. 



(submitted to Physical Review C) 



2 



M. Apollonio/'0 A. Artamonov,^^!! A. Bagulya,^ G. BaiT,^ A. Blondel,^ F. Bobisut^**),^ 
M. Bogomilov/ M. Bonesini,^'EI C. Booth,^ S. Borghi,5'EI S. Bunyatov/o J. Burguet-Castell," 
M.G. Catanesi,i2 A. Cervera-Villanueva," P. Chimenti,i L. Coney,'0 E. Di Capua/^ u. Dore," 
J. Dumarchez,i5 R. Edgecock,!^ M. Ellis/6,0 F. Ferri,* U. Gastaldi,!^ S. Giani,^ G. Giannini,i 
D. Gibin(**),6 S. Gilardoni,^ P. Gorbunov,2.0 C. Gofiling,!^ J.J. Gomez-Cadenas," A. Grant,^ 
J.S. GraulicKiS'S G. Gregoirei^ V. Grichine,^ A. Grossheim,^^!! A. Guglielmi^*),^ L. Hewlett,^ 
A. Ivanchenko,^'! V. Ivanchenko,^'0 A. Kayis-Topaksu,^'@ M. Kirsanov, D. Kolev/ A. Krasnoperov/° 
J. Martm-Albo," C. Meurer,^! M. Mezzetto^*),^ G. B. Mills,'Q M.C. Morone,5'|EI P. Novella," 
D. Orestano^**),^^ V. Palladino,^^ J. Panman,^ I. Papadopoulos,^ F. Pastore(**\^^ S. Piperov,^"^ 
N. Polukliina,3 B. Popov,i"'EI G. Prior,^.^ E. Radicioni,!^ D. Schmitz,'EI R. Schroeter,^ G Skoro,'^ 
M. Sorel," E. Tchemiaev,^ P. Temnikov,^^ V. Tereschenko,i" A. Tonazzo(**\22 l. Tortora(*\22 
R. Tsenov,^ I. Tsukerman,^.^ G. Vidal-Sitjes,i3'0 C. Wiebuscli,2'[p| and P. Zucchelli^'S 

^ Universita degli Studi e Sezione INFN, Trieste, Italy 

^ CERN, Geneva, Switzerland 

^P. N. Lebedev Institute of Physics (FIAN), 
Russian Academy of Sciences, Moscow, Russia 

''Nuclear and Astrophysics Laboratory, University of Oxford, UK 

''Section de Physique, Universite de Geneve, Switzerland 

^Sezione INFN^*'' and Universita degli Studi''**\ Padova, Italy 

'' Faculty of Physics, St. Kliment Ohridski University, 
Sofia, Bulgaria 

^Sezione INFN Milano Bicocca, Milano, Italy 

^ Dept. of Physics, University of Sheffield, UK 

^" Joint Institute for Nuclear Research, JINR Dubna, Russia 

Instituto de Fisica Corpuscular, IFIC, GSIG and Universidad de Valencia, Spain 

''^Sezione INFN, Ban, Italy 

'''^Universita degli Studi e Sezione INFN, Ferrara, Italy 

Universita "La Sapienza" e Sezione INFN Roma I, Roma, Italy 

LPNHE, Universites de Paris VI et VII, Paris, France 

'"Rutherford Appleton Laboratory, Ghilton, Didcot, UK 

Laboratori Nazionali di Legnaro dell' INFN, Legnaro, Italy 

Institut fiir Physik, Universitdt Dortmund, Germany 

"^Institut de Physique Nucleaire, UCL, Louvam-la-Neuve, Belgium 

"^"institute for Nuclear Research, Moscow, Russia 

^'institut fiir Physik, Forschungszentrum Karlsruhe, Germany 

'"^^Sezione INFN^*'' and Universita^**'' Roma Tre, Roma, Italy 

Universita "Federico II" e Sezione INFN, Napoli, Italy 

Institute for Nuclear Research and Nuclear Energy, Academy of Sciences, 

Sofia, Bulgaria 



PACS numbers: 25.40-h, 25.40.Ep 



3 



I. INTRODUCTION 



In many particle and astroparticle physics exper- 
iments the knowledge of liadron production is re- 
quired as an external input to make optimal use 
of the recorded data and to help design the ex- 
perimental facilities. The HARP experiment p] 
is motivated by this need for precise hadron pro- 
duction measurements. It has taken data with 
beams of pions and protons with momenta from 
1.5 GeV/c to 15 GeV/c hitting nuclear targets 
made of a large range of materials. To provide 
a large angular and momentum coverage of the 
produced charged particles the experiment com- 
prises two spectrometers, a forward spectrometer 
built around a dipole magnet covering the angular 
range up to 250 mrad and a large-angle spectrom- 
eter constructed in a solenoidal magnet with an 
angular acceptance of 0.35 rad < 9 < 2.15 rad, 
based on a time-projection-chamber (TPC). 

The main HARP objectives are to measure pion 
yields for a quantitative design of the proton driver 
of future superbeams (high-intensity conventional 
beams) and a neutrino factory Q , to provide mea- 
surements to improve calculations of the atmo- 
spheric neutrino flux [s'-'g'I and to measure particle 
yields as input for the flux calculation of acceler- 
ator neutrino experiments 01 , such as K2K [1, Q , 
MiniBooNE fio"! and SciBooNE [TT|. The momen- 
tum range of the incoming particles presented here 
corresponds to a momentum region of great inter- 
est for neutrino beams and are in a region far from 
coverage by earlier dedicated hadroproduction ex- 
periments [H, [l^. In addition to these specific 
aims, the data provided by HARP are valuable for 
validating hadron production models used in sim- 
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ulation programs. These simulations are playing 
an important role in the interpretation and design 
of modern particle-physics experiments. In par- 
ticular, the simulation of calorimeter response and 
secondary interactions in tracking systems needs to 
be supported by experimental hadron production 
data. Results on forward production of charged pi- 
ons by incident protons are the subject of previous 
HARP publications [Ti-[l8t. The analysis of re- 
sults on charged pion production with charged pion 
beams on the full range of targets can be found in 
Ref. figl ^ 

In this paper, measurements of the double- 
differential cross-section, d'^a^ /dpdil, for forward 
production of protons by incident charged pions or 
protons of 3 GeV/c, 5 GeV/c, 8 GeV/c, 8.9 GeV/c 
(Be only), 12 GeV/c and 12.9 GeV/c (Al only) mo- 
mentum impinging on a thin solid beryllium, car- 
bon, aluminium, copper, tin, tantalum and lead 
targets of 5% nuclear interaction length (Ai) thick- 
ness are presented. 

To our knowledge no high statistics proton pro- 
duction data at low momenta (< 15 GeV/c) in the 
forward direction (< 250 mrad) with incident pro- 
tons or charged pions on nuclear targets have been 
published. 

Data were taken in the T9 beam of the CERN PS. 
The collected statistics, for the different nuclear 
targets, are reported in Table HI 



A. Experimental apparatus 

The HARP experiment makes use of a large- 
acceptance spectrometer consisting of a forward 
and a large-angle detection system. The HARP 
detector is shown in Fig. [1] and is fully described 
in Ref. [13] ■ The forward spectrometer - based on 
five modules of large area drift chambers (NDCl- 
5) and a dipole magnet complemented by a 
set of detectors for particle identification (PID): a 
time-of-fiight wall (TOFW) [il, a large Cherenkov 
detector (CHE) and an electromagnetic calorime- 
ter (ECAL) - covers polar angles up to 250 mrad, 
which is well matched to the angular range of in- 
terest for the measurement of hadron production 
to calculate the properties of conventional neu- 
trino beams. The discrimination power of time-of- 
fiight below 3 GeV/c and the Cherenkov dectector 
above 3 GeV/c are combined to provide power- 



^ results on the large angle charged pion production are 
instead presented in Refs. [20l and [2ll 
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TABLE I: Total number of events and tracks used in the various nuclear 5% Ai target data sets and the number 
of incident protons and charged pious on target as calculated from the pre-scaled incident beam triggers. First 
entries (total DAQ events) are for the positive and negative beam; then numbers are given for incident protons, 
tt"*", 7r~ in units of 10'^ events. 



Data set (GeV/c) 




3 


5 


8 


8.9 


12 


12.9 


beam polarity 




+ 




+ 




+ 




+ 


+ 




+ 


Total DAQ events 


(Be) 


1113 


2233 


1296 


1798 


1935 


1585 


5868 


1207 


1227 






(C) 


1345 


1831 


2628 


1279 


1846 


1399 




1062 


646 






(Al) 


1159 


1523 


1789 


920 


1707 


1059 




619 


741 


4713 




(Cu) 


624 


3325 


2079 


1805 


2089 


1615 




745 


591 






(Sn) 


1637 


1972 


2828 


1625 


2404 


1408 




1803 


937 






(Ta) 


1783 


994 


2084 


1435 


1965 


1505 




866 


961 






(Pb) 


1911 


1282 


2111 


2074 


2266 


1496 




487 


1706 





Data set (GeV/c) 






3 






5 






8 




8 


9 




12 




12 


9 


particle type 




P 


7T + 




P 


7V + 




P 


7T + 


TT 


P 


7r+ 


P 


7r+ 


TV 


P 


7V + 


Acc. beam part. 


(Be) 


99 


246 


731 


289 


384 


914 


761 


341 


826 


2103 


1278 


580 


76 


693 






with forw. int. 


(C) 


101 


257 


299 


542 


754 


530 


709 


358 


772 






470 


41 


352 








(Al) 


86 


213 


486 


376 


523 


308 


637 


335 


611 






306 


27 


435 


2116 


332 




(Cu) 


73 


168 


1185 


408 


611 


850 


741 


397 


966 






363 


33 


347 








(Sn) 


217 


467 


778 


528 


819 


732 


818 


481 


804 






856 


79 


584 








(Ta) 


281 


561 


426 


398 


600 


671 


668 


388 


893 






403 


37 


536 








(Pb) 


310 


611 


473 


387 


594 


997 


758 


444 


896 






221 


20 


839 






Final state p 


(Be) 


7.2 


1.7 


1.5 


18.0 


3.8 


3.1 


37.0 


4.3 


4.8 


86.4 


15.6 


19.3 


1.0 


5.7 






selected with PID 


(C) 


6.2 


1.6 


0.6 


29.1 


6.7 


1.7 


32.2 


4.4 


4.1 






14.8 


0.5 


2.5 








(Al) 


4.2 


1.3 


0.8 


19.0 


4.8 


1.0 


26.8 


3.9 


2.9 






10.6 


0.4 


2.8 


73.3 


5.0 




(Cu) 


2.5 


0.8 


1.6 


18.1 


5.6 


2.9 


31.5 


5.2 


6.0 






12.9 


0.6 


3.1 








(Sn) 


4.9 


1.5 


0.9 


20.8 


6.9 


2.5 


33.6 


6.4 


3.9 






30.2 


1.3 


5.2 








(Ta) 


4.8 


1.4 


0.4 


12.9 


4.6 


2.1 


25.6 


5.0 


4.2 






14.2 


0.6 


4.8 








(Pb) 


4.4 


1.4 


0.3 


12.0 


4.0 


2.8 


26.6 


5.2 


4.4 






7.5 


0.3 


7.4 







ful separation of forward pions and protons [25|. 
The calorimeter is used only for separating pions 
and electrons when characterizing the response of 
the other detectors. The muon contamination of 
the beam is measured with a 1.4 m wide muon 
identifier (BMI). At the downstream end of the 
spectrometer, after a 0.4 m long iron absorber, it 
is made of an iron-scintillator sandwich with five 
planes of six scintillator each, read out at both 
sides, giving a total thickness of 6.4 Aj. 

The large-angle spectrometer consists of a Time 
Projection Chamber (TPC) and Resistive Plate 
Chambers (RPCs), located inside a solenoidal 
magnet. It has a large acceptance in the momen- 
tum and angular range for the pions relevant to 
the production of the muons in a neutrino factory. 
This system is not used in the present analysis. 

The HARP experiment, located in the T9 beam of 
the CERN PS, took data in 2001 and 2002. The 
momentum definition of the T9 beam is known 
with a precision of the order of 1% [26l|. 



The target is placed inside the inner field cage 
(IFC) of the TPC, in an assembly that can be 
moved in and out of the solenoid magnet. The 
solid targets used for the measurements reported 
here have a cylindrical shape with a nominal diam- 
eter of about 30 mm. Their thickness is equivalent 
to about 5% of Ai. 

A set of four multi-wire proportional chambers 
(MWPCs) measures the position and direction of 
the incoming beam particles with an accuracy of 
Ril mm in position and «0.2 mrad in angle per 
projection. A beam time-of-fiight system (BTOF) 
measures the time difference of particles over a 
21.4 m path-length. It is made of two identi- 
cal scintillation hodoscopes, TOFA and TOFB 
(originally built for the NA52 experiment [27]), 
which, together with a small target-defining trigger 
counter (TDS, also used for the trigger), provide 
particle identification at low energies. This pro- 
vides separation of pions, kaons and protons up 
to 5 GeV/c and determines the initial time at the 
interaction vertex (to). The timing resolution of 
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Drift Chambers 




ECAL 



FIG. 1: Schematic layout of the HARP detector. The convention for the coordinate system is shown in the 
lower-right corner. The three most downstream (unlabelled) drift chamber modules are only partly equipped 
with electronics and are not used for tracking. The detector covers a total length of 13.5 m along the beam axis 
and has a maximum width of 6.5 m perpendicular to the beam. The beam muon identifier is visible as the most 
downstream detector (white box). 



the combined BTOF system is about 70 ps. A 
system of two N2-filled Cherenkov detectors (BCA 
and BCB) is used to tag electrons at low ener- 
gies and pions at higher energies. The electron 
and pion tagging efficiency is found to be close to 
100%. At the beam energy used for this analysis 
the Cherenkov counters select all particles lighter 
than protons, while the BTOF is used to reject 
ions. A set of trigger detectors completes the beam 
instrumentation. 

The beam of positive particles used for this mea- 
surement contains mainly positrons, muons, pi- 
ons and protons, with small components of kaons, 
deuterons and heavier ions. Its composition de- 
pends on the selected beam momentum. The pro- 
ton fraction in the incoming positive-particle beam 
varies from - 35% at 3 GeV/c to - 92% at 
12 GeV/c. The negatively-charged particle beam 
is mainly composed of pions with small background 
components of muons and electrons. 

At the first stage of the analysis a favoured beam 
particle type is selected using the beam time-of- 
flight system and the two Cherenkov counters. A 
value of the pulse height consistent with the ab- 
sence of a signal in both beam Cherenkov detec- 
tors distinguishes protons (and kaons) from elec- 
trons and pions. We also ask for time measure- 
ments to be present which are needed for calcu- 
lating the arrival time of the beam proton at the 
target. The beam TOF system is used to reject 
ions, such as deuterons, and to separate protons 



from pions at low momenta. At 3 GeV/c, the 
TOF measurement allows the selection of pions 
from protons to be made at more than Scr. In 
most beam settings the nitrogen pressure in the 
beam Cherenkov counters was too low for kaons to 
be above the threshold. Kaons are thus contained 
in the proton sample. However, the fraction of 
kaons has been measured in the 12.9 GeV/c beam 
configuration and is found to contribute less than 
0.5%, and hence it is negligible both in the pion 
and proton beam sample. Electrons radiate in the 
Cherenkov counters and would be counted as pi- 
ons. In the 3 GeV/c beam, electrons are identi- 
fied by both Cherenkov counters, since the pressure 
was such that pions remained below threshold. In 
the 5 GeV/ c beam electrons could be tagged by one 
Cherenkov counter only, while the other Cherenkov 
counter was used to tag pions. The c/tt fraction 
was measured to be 1% in the 3 GeV/c beam and 
< 10"'^ in the 5 GeV/c beam. By extrapolation 
from the lower-energy beam settings this electron 
contamination can be estimated to be negligible 
(< lO^'^) for the beams where it cannot be mea- 
sured directly. More details on the beain particle 
selection can be found in [T^ and Refs. p2| . 

In addition to the momentum-selected beam of 
protons and pions originating from the T9 produc- 
tion target one expects also the presence of muons 
from pion decay both downstream and upstream 
of the beam momentum selection. Therefore, pre- 
cise absolute knowledge of the muon rate incident 
on the HARP targets is required when measure- 
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ments of particle production with incident pions 
are performed. The particle identification detec- 
tors in the beam do not distinguish muons from 
pions. A separate measurement of the muon com- 
ponent has been performed using data-sets with- 
out target ("empty-target data-sets"). Since the 
empty-target data were taken with the same beam 
parameter settings as the data taken with tar- 
gets, the beam composition can be measured in 
the empty-target runs using the forward spectrom- 
eter and then used as an overall correction for 
the counting of pions in the runs with targets. 
Muons are recognized by their longer range in the 
beam muon identifier (BMI). The punch-through 
background in the BMI is measured counting the 
protons (identified with the beam detectors) thus 
mis-identified as muons by the BMI. A compari- 
son of the punch-through rate between simulated 
incoming pions and protons was used to deter- 
mine a correction for the difference between pi- 
ons and protons and to determine the system- 
atic error. This difference is the dominant sys- 
tematic error in the beam composition measure- 
ment. The aim was to determine the composition 
of the beam as it strikes the target, thus muons 
produced in pion decays after the HARP target 
should be considered as a background to the mea- 
surement of muons in the beam. The rate of these 
latter background muons, which depends mainly 
on the total inelastic cross-section and pion decay, 
was calculated by a Monte Carlo simulation using 
GEANT4 [28]. The muon fraction in the beam 
(at the target) is obtained taking into account the 
efficiency of the BMI selection criteria as well as 
the punch-through and decay backgrounds. The 
analyses for the various beam settings give results 
for R = n/in + tt) of (4.2±1)% and (5.2±1)% for 
the low-momentum beams (3 GeV/c and 5 GeV/c) 
and between (4.1±1)% and (2.8±1)% for the high- 
est momenta (from 8 GeV/c to 12.9 GeV/c). The 
uncertainty in these fractions is dominated by the 
systematic uncertainty in the punch-through back- 
ground. The fact that the background does not 
scale with the decay probability for pions is due to 
the limited acceptance of the beam-line to trans- 
port the decay muons. The muon contamination 
is taken into account in the normalization of the 
pion beam and adds a systematic error of 1% to 
the overall normalization. Only events with a sin- 
gle reconstructed beam track in the four MWPCs, 
good timing measurements in BTOF and no signal 
in the beam halo counters are accepted. 

A downstream trigger in the forward scintillator 
trigger plane (FTP) was required to record an 
event, accepting only tracks with a trajectory out- 
side the central hole (60 mm) which allows beam 
particles to pass. 

The length of the accelerator spill is 400 ms with a 



typical intensity of 15 000 beam particles per spill. 
The average number of events recorded by the data 
acquisition ranges from 300 to 350 per spill. 

The absolute normalization of the cross-section 
was performed using "incident-proton" triggers. 
These are triggers where the same selection on the 
beam particle was applied but no selection on the 
interaction was performed. The rate of this trig- 
ger was down-scaled by a factor 64. These un- 
biased events are used to determine iVpot in the 
cross-section formula (1), see later. 



II. DATA ANALYSIS AND 
CROSS-SECTION CALCULATION. 



This analysis is similar to the one reported in ref- 
erence [131 . For the current analysis we have used 
identical reconstruction and PID algorithms. Sec- 
ondary track selection criteria, described in [T6j . 
are optimized to ensure the quality of momentum 
reconstruction and a clean time-of-flight measure- 
ment while maintaining a high reconstruction effi- 
ciency. 

The background induced by interactions of beam 
particles in the materials outside the target is mea- 
sured by taking data without a target in the target 
holder ("empty target data"). These data are sub- 
ject to the same event and track selection criteria 
as the standard data sets and are subtracted bin- 
by-bin. 

The collected event statistics on the different solid 
targets is summarised in Table HI 

The double-differential cross-section for the pro- 
duction of a particle of type a is calculated as fol- 
lows: 

dpdn^^"' ^' NApt iVpot Ap,Ar2j ^ ' 
where 

• ^p^Q {Pi , dj ) is the cross-section in 
mb/(GeV/c sr) for the particle type a 
( a proton in our case) for each true momen- 
tum and angle bin {pi , 9j ) covered in this 
analysis; 

• (p^, d'j) is the number of particles of type 
a' in bins of reconstructed momentum p'^ and 
angle 9'^ in the raw data, after subtraction 
of empty target data (due to beam protons 
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interacting in material other than the nu- 
clear target). These particles must satisfy 
the event, track and PID selection criteria. 

• J^^iap'e'a' is correction matrix which 
accounts for finite efficiency and resolution 
of the detector. It unfolds the true vari- 
ables Pi ,Oj,a from the reconstructed vari- 
ables p ■ , 9j , a' and corrects the observed 
particle number to take into account ef- 
fects such as reconstruction efficiency, accep- 
tance, absorption, pion decay, tertiary pro- 
duction, PID efficiency and PID misidentifi- 
cation rate. 

• ivl^' iv^ and are normalization 
factors, namely: 

is the number of target nuclei per 
unit area ^; 

A^pot is the number of protons on target 
(particles on target); 

Api and Aflj are the bin sizes in mo- 
mentum and solid angle, respec- 
tively ^. 

We do not make a correction for the attenuation 
of the beam in the target, so that strictly speaking 
the cross-sections are valid for Aj = 5% targets. 

The calculation of the correction matrix 
^PiOjap'e' a' is a rather difficult task. Vari- 
ous techniques are described in the literature to 
obtain this matrix. The method applied here 
and called UFO in fl^ is the unfolding method 
introduced by D'Agostini ^ ^. 

The Monte Carlo simulation of the HARP setup 
is based on GEANT4 [H]. The detector materials 
are accurately described in this simulation as well 



^ A - atomic mass, Nj^ - Avogadro number, p - target den- 
sity and t - target thickness 
3 Api = p™''" - p™'". An J = 27r(cos(ej'"") ~ cos(e™=''')) 
* The unfolding method tries to put in correspondence 
the vector of measured observables (such as particle mo- 
mentum, polar angle and particle type) Xmeas with the 
vector of true values xtruc using a migration matrix: 
aJmoas = A^migr X ^^truc- The goal of the method is to 
compute a transformation (correction matrix) to obtain 
the expected values for xtrue from the measured ones. 
The most simple and obvious solution, based on simple 
matrix inversion M~} , is usually unstable and is domi- 
nated by large variances and strong negative correlations 
between neighbouring bins. In the method of D' Agostini, 
the correction matrix M^^'-' tries to connect the mea- 
surement space (effects) with the space of the true values 
(causes) using an iterative Bayesian approach, based on 
Monte Carlo simulations to estimate the probability for 
a given effect to be produced by a certain cause. 



as the relevant features of the detector response 
and the digitization process. All relevant physics 
processes are considered, including multiple scat- 
tering, energy loss, absorption and re-interactions. 
The track reconstruction used in this analysis and 
the simulation are identical to the ones used for the 
7r+ production in p-Be collisions ^5[]. A detailed 
description of the corrections and their magnitude 
can be found there. 

The reconstruction efficiency (inside the geometri- 
cal acceptance) is larger than 95% above 1.5 GeV/c 
and drops to 80% at 0.5 GeV/c. The requirement 
of a match with a TOFW hit has an efficiency be- 
tween 90% and 95% independent of momentum. 
The electron veto rejects about 1% of the pions 
and protons below 3 GeV/c with a remaining back- 
ground of less than 0.5%. Below Cherenkov thresh- 
old the TOFW separates pions and protons with 
negligible background and an efficiency of «94% 
for protons at low momentum increasing to ~98% 
at threshold. Above Cherenkov threshold the ef- 
ficiency for protons is greater than 98% with less 
than 1% of pions mis-identified as protons. 

The absorption and decay of particles is simulated 
by the Monte Carlo. The absorption correction 
is on average 20%, approximately independent of 
momentum. Uncertainties in the absorption of sec- 
ondaries in the dipole spectrometer material are 
taken into account by a variation of 10% of this 
effect in the simulation. The uncertainty in the 
production of background due to tertiary particles 
is larger. The average correction is ~10% and up to 
20% at 1 GeV/c. The correction includes reinter- 
actions in the detector material as well as a small 
component coming from reinteractions in the tar- 
get. The validity of the generators used in the 
simulation was checked by an analysis of HARP 
data with incoming protons and charged pions on 
aluminium and carbon targets at lower momenta 
(3 GeV/c and 5 GeV/c). A 30% variation of the 
secondary production was applied. The average 
empty-target subtraction amounts to Ri20%. 

Owing to the redundancy of the tracking system 
downstream of the target the detection efficiency 
is very robust under the usual variations of the 
detector performance during the long data tak- 
ing periods. Since the momentum is reconstructed 
without making use of the upstream drift chamber 
module (which is more sensitive in its performance 
to the beam intensity) the reconstruction efficiency 
is uniquely determined by the downstream system. 
No variation of the overall efficiency has been ob- 
served. The performance of the TOFW and CHE 
system have been monitored to be constant for the 
data taking periods used in this analysis. The cal- 
ibration of the detectors was performed on a day- 
by-day basis. 
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A. Error estimation 



The total statistical error of the corrected data is 
composed of the statistical error of the raw data 
and of the statistical error of the unfolding pro- 
cedure, as the unfolding matrix is obtained from 
the data themselves, thus contributing also to the 
statistical error. The statistical error provided by 
the unfolding program is equivalent to the propa- 
gated statistical error of the raw data. In order to 
calculate the statistical error of the unfolding pro- 
cedure a separate analysis is applied, as described 
in [l^,[33l- Its conclusion is that the statistical 
error provided by the unfolding procedure has to 
be multiplied globally by a factor of 2, which is 
done for the analyses described here. This factor 
is somewhat dependent on the shape of the distri- 
butions. 

Different types of sources induce systematic errors 
for the analysis described here: track yield cor- 
rections (~ 5%), particle identification 0.1%), 
momentum and angular reconstruction 1%) ^. 
The strategy to calculate these systematic errors 
and the different methods used for their evalua- 
tion are fully described in [l^ . As a result of these 
systematic error studies, each error source can be 
represented by a covariance matrix. The sum of 
these matrices describes the total systematic error, 
as explained in [l6j . 

The experimental uncertainties are shown for a 
typical target (Be) in Figure [2] at 5 GeV/c and 
12 GeV/c incident beam momenta for incident pro- 
tons and negative pions. They are very similar 
for 7r+ and at the other beam energies. Going 
from lighter (Be, C) to heavier targets (Ta, Pb) 
the corrections for 7r° (conversion) and absorp- 
tion/tertiares increase. 

On average the total integrated systematic error 
is around 5 — 6%, with a differential bin-to-bin 
systematic error of the order of 10 — 11%, to be 
compared with a statistical integrated (bin-to-bin 
differential) error of ~ 2 - 3% (~ 10 - 13%). Sys- 
tematic and statistical errors are roughly of the 
same order. 



^ The quoted error in parenthesis refers to fractional error 
of the integrated cross-section ('5Jjt(%)) in the kinematic 
range covered by the HARP experiment 



III. EXPERIMENTAL RESULTS 



The measured double-differential cross-sections for 
the production of forward protons in the labora- 
tory system as a function of the momentum and 
the polar angle for each incident beam momentum 
are shown in Figures [3] and |4] for two typical solid 
targets: beryllium and lead, as an example of a 
light and a heavy target. The error bars shown 
are the square-roots of the diagonal elements in 
the covariance matrix, where statistical and sys- 
tematic uncertainties are combined in quadrature. 
The correlation of the statistical errors (introduced 
by the unfolding procedure) are typically smaller 
than 20% for adjacent momentum bins and even 
smaller for adjacent angular bins. The correlations 
of the systematic errors are larger, typically 80% 
for adjacent bins. The overall scale error (< 2%) is 
not shown. The results of the analysis for all solid 
targets are fully tabulated in Appendix A ^. 

The dependence of the averaged proton yields 
on the incident beam momentum is shown in 
Fig. [5l The proton yields, averaged over two 
angular regions (0.05 rad < 9 < 0.15 rad and 
0.15 rad < 9 < 0.25 rad) and four momentum re- 
gions (0.5 GeV/c < p < 1.5 GeV/c, 1.5 GeV/c < 
p < 2.5 GeV/c, 2.5 GeV/c < p < 3.5 GeV/c and 
3.5 GeV/c < p < 4.5 GeV/c), are shown. 

The dependence of the averaged proton yields 
on the atomic number A is shown in Fig. [S) 
The yields are shown, averaged over two an- 
gular regions (0.05 rad < 9 < 0.15 rad and 
0.15 rad < 9 < 0.25 rad) and four momentum re- 
gions (0.5 GeV/c < p < 1.5 GeV/c, 1.5 GeV/c < 
p < 2.5 GeV/c, 2.5 GeV/c < p < 3.5 GeV/c and 
3.5 GeV/c < p < 4.5 GeV/c). 

The particle production ratios /tt^ and tt^/p, 
in the integrated angular range 0.05 rad < 9 < 
0.25 rad, are reported instead in Figures [7] and [5] 
for two typical targets: beryllium and tantalum 
and the different incident beam particles. They 
clearly reflect the charge polarity of the incident 
beam particle, especially at high secondary mo- 
menta. 



the scale error has not been included to make it possible 
to calculate e.g. integrated particle ratios taking it into 
account only when applicable, i.e. when different beams 
are compared. 
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A. Comparison with Monte Carlo generators. 

In the following we will show only some compar- 
isons with two widely used Monte Carlo simulation 
packages: MARS ^ and GEANT4 H ^ using 
different generator models. The comparison will 
be shown for a limited set of plots and only for the 
Be and Ta targets, as examples of a light and a 
heavy target. In both generators, no single model 
is applicable to all energies and a transition be- 
tween low energy models and high energy models, 
at about 5-10 GeV, is needed. 

The lack of hadron nuclei collisions data with small 
errors (both statistics and systematics) on an ex- 
tended set of thin targets has been, up to now, an 
obstacle for a serious tuning of these models. Ded- 
icated simulations, such as GiBUU [32], may also 
profit from the availability of our data. 

At intermediate energies (up to 5-10 GeV), 
GEANT4 uses two types of intra-nuclear cascade 
models: the Bertini model [13, 01 (valid up to 
~ 10 GeV) and the binary model [35| (valid up to 
~ 3 GeV). Both models treat the target nucleus in 
detail, taking into account density variations and 
tracking in the nuclear field. The binary model is 
based on hadron collisions with nucleons, giving 
resonances that decay according to their quantum 
numbers. The Bertini model is based on the cas- 
cade code reported in [36l | and hadron collisions are 
assumed to proceed according to free-space partial 
cross sections and final state distributions mea- 
sured for the incident particle types. Details of 
the nuclear density and the Pauli blocking are then 
taken into account. 

At higher energies, instead, two parton stringmod- 
els, the quark-gluon string (QGS) model [H, \3l\ 
and the Fritiof (FTP) model [33l are used, in addi- 
tion to a High Energy Parametrized model (HEP) 
derived from the high energy part of the GHEISHA 
code used inside GEANT3 |38i| . The parametrized 
models of GEANT4 (HEP and LEP) are intended 
to be fast, but conserve energy and momentum on 
average and not event-by-event. 

A realistic GEANT4 simulation is built by com- 
bining models and physics processes into what is 
called a "physics list" , that is included in the stan- 
dard GEANT4 Toolkit release. Each physics list 
corresponds to a collection of models suitable for 
a given user problem. 

As examples, the QGSP physics list is based on 



the QGS model, the pre-compound nucleus model 
and some of the Low Energy Parametrized (LEP) 
model while the LHEP physics list [39| is based 
on the parametrized LEP model and HEP mod- 
els. Currently the most widely used physics list 
in LHC experiments is the so-called QGSP-Bert 
physics list, see reference for details. 

The MARS code system [sij uses as basic model an 
inclusive approach multi-particle production orig- 
inated by R. Feynman. Above 5 GeV phenomeno- 
logical particle production models are used, while 
below 5 GeV a cascade-exciton model [41| com- 
bined with the Fermi break-up model, the coales- 
cence model, an evaporation model and a multi- 
fragmentation extension are used instead. 

The comparison, just outlined in our paper, be- 
tween data and models is reasonable, but some dis- 
crepancies are evident for some models especially 
at lower energies and small angles. Discrepancies 
up to a factor of three are seen. 

The full set of HARP data, taken with targets 
spanning the full periodic table of elements, with 
small total errors and full coverage of the solid an- 
gle in a single detector may help the validation of 
models used in hadronic simulations in the diffi- 
cult energy range between 3 and 15 GeV/c of inci- 
dent momentum, as done e.g. in reference (32j for 
charged pion production. 



The GEANT4 version used is 9.2p01 



Also this model, at low energy, has its root in the 
GHEISHA code inside GEANT3. 
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FIG. 9: Comparison of HARP double-differential pro- 
ton cross sections for p-Be, 7r~-Be, 7r"'"-Be interactions 

at 3 GcV/c with GEANT4 and MARS MC predictions, 
using several generator models (see text for details). 



FIG. 10: Comparison of HARP double-differential 
proton cross sections for p-Be, 7r~-Be, 7r"'"-Be inter- 
actions at 5 GeV/c with GEANT4 and MARS MC 
predictions, using several generator models (see text 
for details). 
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FIG. 11: Comparison of HARP double-differential 
proton cross sections for p-Be, 7r~-Be, 7r"'"-Be inter- 
actions at 8 GeV/c with GEANT4 and MARS MC 
predictions, using several generator models (see text 
for details). 



FIG. 12: Comparison of HARP double-differential 
proton cross sections for p-Be, 7r~-Be, 7r"'"-Be inter- 
actions at 12 GcV/c with GEANT4 and MARS MC 
predictions, using several generator models (see text 
for details). 
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FIG. 13: Comparison of HARP double-differential 
proton cross sections for p-Ta, 7r~-Ta, ir'^-Ta, inter- 
actions at 3 GcV/c with GEANT4 and MARS MC 
predictions, using several generator models (see text 
for details). 



FIG. 14: Comparison of HARP double-differential 
proton cross sections for p-Ta, 7r~-Ta, 7r"'"-Ta inter- 
actions at 5 GeV/c with GEANT4 and MARS MC 
predictions, using several generator models (see text 
for details). 
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FIG. 15: Comparison of HARP double-differential 
proton cross sections for p-Ta, 7r~-Ta, ir'^-Ta, inter- 
actions at 8 GcV/c with GEANT4 and MARS MC 
predictions, using several generator models (see text 
for details). 



FIG. 16: Comparison of HARP double-differential 
proton cross sections for p-Ta, 7r~-Ta, 7r"'"-Ta inter- 
actions at 12 GeV/c with GEANT4 and MARS MC 
predictions, using several generator models (see text 
for details). 
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IV. SUMMARY AND CONCLUSIONS 



In this paper we report our results on double- 
differential cross sections for the forward produc- 
tion of protons in the kinematic range 0.5 GeV/c< 
Ptt < S> GeV/c and 0.05 rad < 6^^ < 0.25 rad from 
the collisions of protons and charged pions of 3, 5, 
8 and 12 GeV/c on beryllium, carbon, aluminium, 
copper, tin, tantalum and lead targets of 5% Ai 
thickness. 

The proton yield averaged over different momen- 
tum and angular ranges increases smoothly with 
the atomic number A of the target and with the 
energy of the incoming beam. 

Comparisons with GEANT4 and MARS genera- 
tors are presented. 

We stress that the HARP data presented here are 
the first precision measurements of forward pro- 
ton production in this kinematic region and may 
have a major impact on the tuning of Monte Carlo 
generators. 
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Appendix A: Cross-section data 

The following tables report the measured differ- 
ential cross-section for forward proton production 
in interactions of 3, 5, 8 and 12 GeV/c momen- 
tum charged pions or protons on different types 
of nuclear targets. The data arc presented in the 
kinematic range of 0.5 GeV/c < < 8 GeV/c 
and 0.05 rad < < 0.25 rad. The overall normal- 
ization uncertainty (< 2%) is not included in the 
reported errors. 



TABLE II: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, dFa^/{dpdQ), for 7r~-Be interactions at 
3,5,8,12 GeV/c. Each row refers to a different 
(Pmin < -P < Pmax,^'rnin < 9 < ^'max) bin, where p 
and 9 are the outgoing proton momentum and polar 
angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance 
matrix are given. 
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TABLE III: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, (fa^/{dpdQ), for Tr^-Be interactions at 
3,5,8,8.9 GeV/c. Each row refers to a different 
(Pmin < -P < Pmax,^rnin < 9 < ^'max) bin, where p 
and 9 are the outgoing proton momentum and polar 
angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance 
matrix are given. 
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TABLE IV: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, (f a'' / {dpdfl) , for p-Be interactions at 
3,5,8,8.9,12 GeV/c. Each row refers to a differ- 
ent (Pjnin < P < Pmax,6'niin < ^ < ^max) bin, where 
p and 6 are the outgoing proton momentum and polar 
angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance 
matrix are given. 
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TABLE V: HARP results for the double-differential 
p production cross-section in the laboratory system, 
dPa"' /{dpdO), for tx~-G interactions at 3,5,8,12 GeV/c. 
Each row refers to a different (p^ain ^ P < 
Pmaxj^'inin ^ S < ^max) bin, whore p and 6 arc the 
outgoing proton momentum and polar angle, respec- 
tively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 



e ■ 

mill 
(rad) 


Omax Pmin Pmax 
(rad) (GeV/c) (GeV/c) 


3 GeV/c 


tfa^/idpdn) 

(barn/(sr GeV/c)) 
5 GeV/c 8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.11±0.02 


0.08±0.02 


0.08±0.02 


0.05±0.02 






1.00 


1.50 


0.056±0.014 


0.068±0.011 


0.057±0.012 


0.060±0.013 






1.50 


2.00 


0.054±0.011 


0.069±0.010 


0.077±0.012 


0.05±0.03 






2.00 


2.50 


0.031±0.008 


0.045±0.007 


0.051±0.010 


0.005±0.006 






2.50 


3.00 


± 


0.026±0.005 


0.052±0.007 


0.035±0.011 






3.00 


3.50 


± 


0.018±0.005 


0.041±0.007 


0.056±0.009 






3.50 


4.00 


± 


0.012±0.003 


0.043±0.006 


0.052±0.007 






4.00 


5.00 


± 


0.008±0.003 


0.020±0.004 


0.055±0.006 






5.00 


6.50 


± 


± 


0.007±0.002 


0.026±0.004 






6.50 


8.00 


± 


± 


± 


0.011±0.002 


0.100 


0.150 


0.50 


1.00 


0.08±0.02 


0.07±0.02 


0.07±0.02 


0.05±0.02 






1.00 


1.50 


0.07±0.02 


0.077±0.013 


0.069±0.012 


0.065±0.014 






1.50 


2.00 


0.071±0.015 


0.060±0.010 


0.061±0.011 


0.036±0.011 






2.00 


2.50 


0.022±0.007 


0.047±0.008 


0.048±0.008 


0.034±0.012 






2.50 


3.00 


± 


0.029±0.006 


0.051±0.006 


0.036±0.008 






3.00 


3.50 


± 


0.011±0.003 


0.025±0.006 


0.048±0.008 






3.50 


4.00 


± 


0.011±0.003 


0.022±0.004 


0.047±0.007 






4.00 


5.00 


± 


0.005±0.002 


0.011±0.003 


0.031±0.005 






5.00 


6.50 


± 


± 


0.003±0.001 


0.009±0.002 






6.50 


8.00 


± 


± 


± 


0.003±0.001 


0.150 


0.200 


0.50 


1.00 


0.09±0.03 


0.08±0.02 


0.06±0.02 


0.07±0.02 






1.00 


1.50 


0.036±0.010 


0.048±0.009 


0.044±0.009 


0.041±0.010 






1.50 


2.00 


0.037±0.010 


0.060±0.010 


0.053±0.008 


0.047±0.012 






2.00 


2.50 


0.031±0.009 


0.034±0.008 


0.048±0.010 


0.022±0.008 






2.50 


3.00 


± 


0.021±0.005 


0.025±0.005 


0.033±0.008 






3.00 


3.50 


± 
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0.021±0.005 


0.041±0.009 






3.50 
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0.013±0.003 


0.019±0.006 






4.00 


5.00 


± 


± 
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0.007±0.002 






5.00 


6.50 


± 


± 


± 


0.002±0.001 






6.50 


8.00 


± 


± 


± 


± 


0.200 


0.250 


0.50 


1.00 


0.09±0.03 


0.09±0.02 


0.06±0.02 


0.05±0.02 






1.00 


1.50 


0.11±0.03 


0.07±0.02 


0.09±0.02 


0.06±0.02 






1.50 


2.00 


0.05±0.02 


0.064±0.014 


0.06±0.02 


0.05±0.02 






2.00 


2.50 


0.034±0.013 


0.036±0.009 


0.08±0.02 


0.04±0.02 






2.50 


3.00 


± 


0.022±0.007 


0.035±0.008 


0.030±0.010 






3.00 


3.50 


± 


0.009±0.004 


0.013±0.004 


0.015±0.006 






3.50 


4.00 
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0.002±0.001 


0.010±0.003 


0.007±0.003 






4.00 


5.00 


± 
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0.005±0.002 


0.004±0.002 






5.00 


6.50 


± 
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O.OOliO.OOl 


± 






6.50 


8.00 


± 


± 


± 


± 
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TABLE VI: HARP results for the double-differential 
p production cross-section in the laboratory system, 
<fa^/{dpdQ), for 7r+-C interactions at 3,5,8,12 GeV/c. 
Each row refers to a different (Pj^i^ < P < 
Pmaxj^'inin ^ S < ^max) bin, whore p and 6 arc the 
outgoing proton momentum and polar angle, respec- 
tively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 



e ■ 

mill 
(rad) 


Omax Pmin Pmax 
(rad) (GeV/c) (GeV/c) 


3 GeV/c 


tfa^/idpdn) 

(barn/(sr GeV/c)) 
5 GeV/c 8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.10±0.04 


0.07±0.02 


0.06±0.02 


0.06±0.05 






1.00 


1.50 


0.09±0.03 


0.08±0.02 


0.07±0.02 


0.07±0.04 






1.50 


2.00 


± 


0.07±0.02 


0.08±0.02 


0.01±0.03 






2.00 


2.50 


± 


0.015±0.016 


0.063±0.014 


0.03±0.08 






2.50 


3.00 


± 


± 


0.036±0.013 


0.055±0.084 






3.00 


3.50 


± 


± 


0.002±0.008 


0.026±0.039 






3.50 


4.00 


± 


± 


0.05±0.02 


0.001±0.012 






4.00 


5.00 


± 


± 


± 


± 






5.00 


6.50 


± 


± 


± 


± 






6.50 


8.00 


± 


± 


± 


0.013±0.008 


0.100 


0.150 


0.50 


1.00 


0.09±0.04 


0.07±0.02 


0.06±0.02 


0.10±0.07 






1.00 


1.50 


0.10±0.03 


0.06±0.01 


0.06±0.02 


± 






1.50 


2.00 


0.05±0.03 


0.06±0.02 


0.06±0.02 


0.09±0.05 






2.00 


2.50 


± 


0.03±0.02 


0.039±0.012 


0.03±0.03 






2.50 


3.00 


± 


± 


0.040±0.010 


0.01±0.02 






3.00 


3.50 


± 


± 


0.041±0.011 


0.024±0.024 






3.50 


4.00 


± 


± 


0.018±0.007 


± 






4.00 


5.00 


± 


± 


0.011±0.005 


± 






5.00 


6.50 


± 


± 


0.006±0.019 


± 






6.50 


8.00 


± 


± 


± 


± 


0.150 


0.200 


0.50 


1.00 


0.10±0.04 


0.10±0.02 


0.05±0.02 


0.09±0.07 






1.00 


1.50 


0.08±0.02 


0.070±0.013 


0.051±0.014 


0.02±0.03 






1.50 


2.00 


0.018±0.023 


0.058±0.011 


0.036±0.012 


0.05±0.04 






2.00 


2.50 


± 


0.029±0.011 


0.029±0.010 


0.018±0.024 






2.50 


3.00 


± 


0.012±0.008 


0.020±0.007 


0.014±0.022 






3.00 


3.50 


± 


0.005±0.006 


0.011±0.005 


0.020±0.046 






3.50 


4.00 


± 


0.011±0.007 


0.015±0.006 


0.005±0.022 






4.00 


5.00 


± 


0.006±0.006 


0.006±0.003 


0.005±0.022 






5.00 


6.50 


± 


± 


O.OOliO.OOl 


0.001±0.006 






6.50 


8.00 


± 


± 


± 


± 


0.200 


0.250 


0.50 


1.00 


0.09±0.04 


0.07±0.03 


0.05±0.02 


0.06±0.06 






1.00 


1.50 


0.04±0.02 


0.06±0.02 


0.09±0.03 


0.03±0.06 






1.50 


2.00 


0.04±0.02 


0.036±0.012 


0.07±0.02 


0.02±0.05 






2.00 


2.50 


0.02±0.02 


0.039±0.013 


0.028±0.013 


0.06±0.09 






2.50 


3.00 


± 


0.024±0.007 


0.019±0.008 


0.01±0.02 






3.00 


3.50 


± 


0.012±0.005 


0.024±0.008 


0.029±0.058 






3.50 


4.00 


± 


0.007±0.003 


0.015±0.005 


0.012±0.025 






4.00 


5.00 


± 


0.004±0.002 


0.014±0.006 


0.002±0.011 






5.00 


6.50 


± 


± 


0.004±0.003 


0.000±0.004 






6.50 


8.00 


± 


± 


± 


0.000±0.005 



TABLE VII: HARP results for the double-differential 
p production cross-section in the laboratory system, 
dPa"' /{dpdO), for p-G interactions at 3,5,8,12 GeV/c. 
Each row refers to a different (p^ain ^ P < 
Pmaxj^'inin ^ S < ^max) bin, whore p and 6 arc the 
outgoing proton momentum and polar angle, respec- 
tively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 



e ■ 

mill 


Omax 


Pmin 


pmax 








(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.12±0.05 


1 3-1-0 03 


10+0 02 


0.09+0.02 






1.00 


1.50 


0.21±0.04 


n 1 3+n 09 

U . -LO 1 U.UZi 


1 99+0 01 5 

yj . ± yj.yjA.yj 


0.088+0.015 






1.50 


2.00 


0.31±0.05 


16±0 02 


16+0 02 


0.090+0.015 






2.00 


2.50 


0.27±0.09 


22+0 02 


16+0 02 


0.12+0.02 






2.50 


3.00 


± 


33±0 03 


20+0 02 


0.15+0.01 






3.00 


3.50 


± 


41+0 03 

W .t: J \—\J . L/t_/ 


27+0 02 


0.14+0.02 






3.50 


4.00 


± 


40+0 03 


30+0 02 


0.18+0.02 






4.00 


5.00 


± 


n 50+0 07 


33+0 02 


0.196+0.013 






5.00 


6.50 


± 


-I- 


24+0 01 

KJ . _1_ KJ .\J ±. 


0.191+0.010 






6.50 


8.00 


± 


-I- 


+ 


0.159+0.008 


0.100 


0.150 


0.50 


1.00 


0.12±0.05 


1 3+0 03 


1 3+0 02 


0.13+0.03 






1.00 


1.50 


0.22±0.04 


19+0 09 


11+0 09 


0.14+0.02 






1.50 


2.00 


0.32±0.05 


n 18+n 02 


1 2+0 02 

yj w A. ^ 1 yj *yj ^ 


0.088+0.014 






2.00 


2.50 


0.30±0.05 


lQ+0 02 


14+0 02 

yj . yj^ 


0.080+0.013 






2.50 


3.00 


± 


22+0 02 


18+0 02 


0.106+0.014 






3.00 


3.50 


± 


22+0 02 


177+0 015 

yj . A- 1 1 1 yj.yjA.ij 


0.099+0.012 






3.50 


4.00 


± 


18+0 02 


155+0 012 

\J- A.yJyJ—l—\J.yj J- 


0.117+0.012 






4.00 


5.00 


± 


1 Q+0 03 


1 37+0 01 1 


0.112+0.010 






5.00 


6.50 


± 


-l- 


068+0 006 


0.064+0.006 






6.50 


8.00 


± 






0.038+0.004 


0.150 


0.200 


0.50 


1.00 


0.08±0.04 


1 3+0 03 


07+0 02 

yj.KJi .s^Kj.yj^ 


0.10+0.02 






1.00 


1.50 


0.17±0.04 


1 Q-l-0 09 


085+0 01 3 


0.063+0.012 






1.50 


2.00 


0.15±0.03 


n 1 4+n 09 


088+0 01 9 


0.076+0.014 






2.00 


2.50 


0.11±0.03 


0.14±0.02 


0.111+0.013 


0.061+0.011 






2.50 


3.00 


± 


0.122±0.014 


0.098+0.011 


0.063+0.010 






3.00 


3.50 


± 


0.093±0.012 


0.083+0.009 


0.054+0.008 






3.50 


4.00 


± 


0.079+0.011 


0.067+0.007 


0.042+0.006 






4.00 


5.00 


± 


0.058+0.007 


0.052+0.006 


0.034+0.005 






5.00 


6.50 


± 


± 


0.015+0.003 


0.015+0.003 






6.50 


8.00 


± 


± 


+ 


0.004+0.001 


0.200 


0.250 


0.50 


1.00 


0.15±0.07 


0.09+0.03 


0.12+0.03 


0.10+0.03 






1.00 


1.50 


0.14±0.04 


0.10+0.02 


0.10+0.02 


0.07+0.02 






1.50 


2.00 


0.09±0.03 


0.07+0.02 


0.07+0.02 


0.036+0.011 






2.00 


2.50 


0.08±0.03 


0.059+0.013 


0.09+0.02 


0.054+0.015 






2.50 


3.00 


± 


0.056+0.011 


0.079+0.014 


0.019+0.006 






3.00 


3.50 


± 


0.031+0.007 


0.060+0.009 


0.020+0.005 






3.50 


4.00 


± 


0.026+0.005 


0.043+0.008 


0.022+0.005 






4.00 


5.00 


± 


0.026+0.005 


0.027+0.005 


0.014+0.004 






5.00 


6.50 


± 


± 


0.007+0.002 


0.005+0.002 






6.50 


8.00 


± 


± 


+ 


0.002+0.001 
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TABLE VIII: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, (f(j^/{dpdQ), for 7r~-Al interactions at 
3,5,8,12 GeV/c. Each row refers to a different 
(Pmin < -P < Pmax,^'rnin < 9 < ^'max) bin, where p 
and 9 are the outgoing proton momentum and polar 
angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance 
matrix are given. 



9 ■ 
mill 

(rad) 


0max Prnin Pmax 
(rad) (GeV/c) (GeV/c) 


3 GeV/c 


<fa'^/{dp<m) 
(barn/(sr GeV/c)) 
5 GeV/c 8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.16±0.03 


0.15±0.04 


0.15±0.03 


0.12±0.03 






1.00 


1.50 


0.12±0.02 


0.13±0.02 


0.14±0.02 


0.11±0.02 






1.50 


2.00 


0.075±0.013 


0.12±0.02 


0.14±0.02 


0.08±0.02 






2.00 


2.50 


0.043±0.010 


0.09±0.02 


0.08±0.02 


0.010±0.013 






2.50 


3.00 


± 


0.047±0.011 


0.104±0.013 


0.08±0.02 






3.00 


3.50 


± 


0.016±0.006 


0.067±0.011 


0.096±0.013 






3.50 


4.00 


± 


0.033±0.009 


0.067±0.009 


0.086±0.011 






4.00 


5.00 


± 


0.012±0.006 


0.045±0.007 


0.104±0.010 






5.00 


6.50 


± 


± 


0.010±0.003 


0.042±0.006 






6.50 


8.00 


± 


± 


± 


0.022±0.004 


0.100 


0.150 


0.50 


1.00 


0.20±0.04 


0.18±0.04 


0.17±0.03 


0.14±0.03 






1.00 


1.50 


0.14±0.02 


0.19±0.03 


0.15±0.02 


0.15±0.02 






1.50 


2.00 


0.09±0.02 


0.16±0.03 


0.12±0.02 


0.09±0.02 






2.00 


2.50 


0.047±0.010 


0.10±0.02 


0.11±0.02 


0.06±0.02 






2.50 


3.00 


± 


0.048±0.011 


0.087±0.011 


0.070±0.013 






3.00 


3.50 


± 


0.021±0.007 


0.056±0.011 


0.097±0.015 






3.50 


4.00 


± 


0.017±0.006 


0.039±0.007 


0.082±0.012 






4.00 


5.00 


± 


0.008±0.004 


0.025±0.005 


0.053±0.008 






5.00 


6.50 


± 


± 


0.005±0.002 


0.015±0.003 






6.50 


8.00 


± 


± 


± 


0.004±0.001 


0.150 


0.200 


0.50 


1.00 


0.16±0.04 


0.09±0.03 


0.16±0.03 


0.15±0.03 






1.00 


1.50 


0.08±0.02 


0.12±0.02 


0.10±0.02 


0.08±0.02 






1.50 


2.00 


0.061±0.013 


0.10±0.02 


0.12±0.02 


0.09±0.02 






2.00 


2.50 


0.052±0.011 


0.08±0.02 


0.08±0.02 


0.07±0.02 






2.50 


3.00 


± 


0.027±0.008 


0.068±0.012 


0.054±0.011 






3.00 


3.50 


± 


0.029±0.009 


0.044±0.010 


0.049±0.009 






3.50 


4.00 


± 


0.002±0.002 


0.026±0.006 


0.045±0.008 






4.00 


5.00 


± 


O.OOOiO.OOl 


0.011±0.003 


0.024±0.005 






5.00 


6.50 


± 


± 


O.OOliO.OOl 


0.004±0.001 






6.50 


8.00 


± 


± 


± 


± 


0.200 


0.250 


0.50 


1.00 


0.16±0.04 


0.16±0.05 


0.18±0.04 


0.17±0.04 






1.00 


1.50 


0.14±0.03 


0.18±0.04 


0.21±0.04 


0.15±0.03 






1.50 


2.00 


0.09±0.02 


0.13±0.03 


0.16±0.04 


0.11±0.03 






2.00 


2.50 


0.030±0.012 


0.09±0.02 


0.13±0.03 


0.08±0.03 






2.50 


3.00 


± 


0.06±0.02 


0.06±0.02 


0.06±0.02 






3.00 


3.50 


± 


0.019±0.009 


0.024±0.008 


0.036±0.009 






3.50 


4.00 


± 


0.005±0.003 


0.008±0.004 


0.021±0.006 






4.00 


5.00 


± 


0.001±0.002 


0.003±0.002 


0.012±0.005 






5.00 


6.50 


± 


± 


0.000±0.002 


O.OOliO.OOl 






6.50 


8.00 


± 


± 


± 


O.OOOiO.OOl 



TABLE IX: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, (f a'' / {dpdO,) , for tt^-AI interactions at 
3,5,8,12,12.9 GeV/c. Each row refers to a differ- 
ent (Pjnin < P < Pmax,^'niin < ^ < ^max) bin, where 
p and 9 are the ougoing proton rnorncnturn and polar 
angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance 
matrix are given. 



e ■ 

mill 


Omax 


Pmin 


pmax 






dPa'^/idpdO) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


12.9 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.16±0.06 


f) 1 7-l-f) QA 


f) 14+0 04 


n 17+0 1Q 


0.18+0.05 






1.00 


1.50 


0.16±0.05 


1 7+0 OS 


1 3+0 03 


10+0 09 


0.12+0.03 






1.50 


2.00 


± 


13±0 03 


11+0 03 


05+0 08 

yj *yj*-j —I — yj *yjyj 


0.14+0.03 






2.00 


2.50 


± 


OQ+n 03 


oq+0 02 


07+0 12 

\j .yj 1 —1— \j . 


0.13+0.02 






2.50 


3.00 


± 


± 


07+0 02 


07+0 08 


0.12+0.02 






3.00 


3.50 


± 


-I- 


03+0 02 


08+0 10 

\j . yj'jf 1 yj . J- yj 


0.11+0.02 






3.50 


4.00 


± 


-I- 


+ 


02+0 08 


0.10+0.02 






4.00 


5.00 


± 


-I- 


+ 


f) 03+0 1 5 


0.14+0.02 






5.00 


6.50 


± 


-l- 




005+0 099 
u.yjyjO-i—yj.yjty.^ 


0.075+0.010 






6.50 


8.00 


± 


-I- 


+ 


n 002+0 127 


0.039+0.007 


0.100 


0.150 


0.50 


1.00 


0.21±0.07 


21-1-0 05 


n 24+ n 06 


01 +0 03 


0.14+0.05 






1.00 


1.50 


0.15±0.04 


1 3-1-0 03 


1 3+0 03 


1 1 +0 08 

W . J- J — 1_ yJ • \JyJ 


0.14+0.03 






1.50 


2.00 


0.07±0.05 


n 1 3+n 03 


10+0 03 


07+0 08 


0.11+0.03 






2.00 


2.50 


± 


05+0 07 


08+0 02 


06+0 06 

\j • yj\j.j—\j • yjyj 


0.09+0.02 






2.50 


3.00 


± 


± 


09+0 02 


07+0 07 

KJ myj t KJ myj t 


0.08+0.02 






3.00 


3.50 


± 




05+0 02 


02+0 03 

\j . yj \ yj . yjrj 


0.09+0.02 






3.50 


4.00 


± 


-I- 


019+0 009 


n 0S9+0 045 


0.075+0.013 






4.00 


5.00 


± 


-I- 


01 7+0 007 


n 064+0 075 


0.068+0.012 






5.00 


6.50 


± 


-I- 


01 1 +0 OOfi 


01 1 +0 026 


0.030+0.006 






6.50 


8.00 


± 


-|- 




005+0 016 

yj •yjyjyj i yj .yj 


0.010+0.003 


0.150 


0.200 


0.50 


1.00 


0.19±0.07 


n 20+n 04 


20+0 05 


10+0 14 


0.08+0.04 






1.00 


1.50 


0.18±0.04 


n io-|-n 09 


19+0 03 


f) iK-i-n 1 1 

\J. J-U^tW. J- J- 


0.12+0.03 






1.50 


2.00 


0.05±0.04 


n 1 1 +n 09 

. -L J \—\J .\J ^ 


07+0 02 


1 5+0 1 1 


0.08+0.02 






2.00 


2.50 


± 


0.08±0.02 


0.06+0.02 


0.06+0.07 


0.07+0.02 






2.50 


3.00 


± 


0.02±0.02 


0.08+0.02 


0.013+0.032 


0.051+0.015 






3.00 


3.50 


± 


0.001±0.006 


0.035+0.012 


0.051+0.076 


0.049+0.014 






3.50 


4.00 


± 


0.013+0.011 


0.013+0.008 


0.028+0.044 


0.037+0.010 






4.00 


5.00 


± 


0.012+0.013 


0.009+0.008 


0.011+0.024 


0.018+0.006 






5.00 


6.50 


± 


± 


0.002+0.002 


0.006+0.015 


0.010+0.003 






6.50 


8.00 


± 


+ 


+ 


+ 


0.002+0.001 


0.200 


0.250 


0.50 


1.00 


0.13±0.07 


0.17+0.06 


0.17+0.05 


0.10+0.17 


0.17+0.06 






1.00 


1.50 


0.14±0.05 


0.11+0.03 


0.11+0.03 


0.08+0.13 


0.08+0.02 






1.50 


2.00 


0.06±0.03 


0.06+0.02 


0.10+0.03 


0.02+0.09 


0.06+0.03 






2.00 


2.50 


0.04±0.03 


0.06+0.02 


0.09+0.03 


0.13+0.22 


0.06+0.02 






2.50 


3.00 


± 


0.042+0.014 


0.06+0.02 


0.01+0.02 


0.04+0.02 






3.00 


3.50 


± 


0.019+0.009 


0.04+0.02 


0.079+0.116 


0.031+0.009 






3.50 


4.00 


± 


0.012+0.006 


0.026+0.009 


0.004+0.016 


0.019+0.008 






4.00 


5.00 


± 


0.009+0.004 


0.021+0.009 


0.000+0.004 


0.008+0.005 






5.00 


6.50 


± 


± 


0.008+0.008 


0.001+0.013 


0.002+0.002 






6.50 


8.00 


± 


± 


± 


0.000+0.004 


0.001+0.001 



TABLE X: HARP results for the double-differential 
p production cross-section in the laboratory 
system, (fi / {dpdKl) , for p-Al interactions at 
3,5,8,12,12.9 GeV/c. Each row refers to a different 
(Pmin < -P < Pmax,^'rnin < 9 < ^'max) bin, where p 
and 9 are the outgoing proton momentum and polar 
angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance 
matrix are given. 



9 ■ 
mill 


0max 


Pmin 


pmax 






(fa^'/idtpdn) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


12.9 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.20±0.08 


f) 22-1-0 05 


f) 22+0 04 


22+0 05 


0.23+0.03 






1.00 


1.50 


0.32±0.07 


n 27+0 03 


24+0 03 


23+0 04 


0.21+0.02 






1.50 


2.00 


0.41±0.08 


33±0 04 


0.23+0.02 


26+0 04 


0.21+0.02 






2.00 


2.50 


0.39±0.12 


n 34+0 04 


26+0 03 


1 7+0 03 


0.22+0.02 






2.50 


3.00 


± 


53±0 04 


35+0 03 

\J t^KJ ^^\J rnXJiJ 


26+0 03 


0.25+0.02 






3.00 


3.50 


± 


n 53+n 05 


37+0 02 


28+0 03 


0.26+0.02 






3.50 


4.00 


± 


n 59+0 05 


35+0 02 


30+0 03 


0.30+0.02 






4.00 


5.00 


± 


n fi2+n OS 


n 41 +n 02 


34+0 03 


0.34+0.02 






5.00 


6.50 


± 


-I- 


n 35+0 02 

\J 'fJ^J —l—\J -yJ^ 


31+0 02 

\J .rj J \_\J . WZj 


0.312+0.011 






6.50 


8.00 


± 


-I- 


+ 


29+0 02 


0.272+0.010 


0.100 


0.150 


0.50 


1.00 


0.28±0.09 


26+0 05 


2.3+0 04 


24+0 05 


0.24+0.03 






1.00 


1.50 


0.30±0.07 




18+0 03 


2fi+0 04 


0.22+0.02 






1.50 


2.00 


0.53±0.09 


n 27+n 03 


24+0 03 

yj . [ \j .yjrj 


22+0 03 

\j * 1 \j *yjfj 


0.21+0.02 






2.00 


2.50 


0.44±0.08 


35+0 04 

\j .tj<j .yj^ 


23+0 03 


18+0 03 


0.21+0.02 






2.50 


3.00 


± 




29+0 03 


20+0 03 


0.20+0.02 






3.00 


3.50 


± 


n 34+n 04 


28+0 02 


20+0 03 


0.20+0.02 






3.50 


4.00 


± 


32+0 03 


22+0 02 


1 7+0 02 


0.189+0.015 






4.00 


5.00 


± 


n 32+0 05 


20+0 02 


20+0 02 


0.190+0.013 






5.00 


6.50 


± 


-I- 


116+0 009 


1 1 7+0 01 2 


0.121+0.009 






6.50 


8.00 


± 


-|- 


+ 


076+0 009 

yj -yj 1 yj —i— yj •yjyj 


0.063+0.006 


0.150 


0.200 


0.50 


1.00 


0.19±0.08 


27+0 06 

I .J.^\J .\J\J 


1 8+0 04 


22+0 05 


0.21+0.03 






1.00 


1.50 


0.27±0.06 


n 21 +n 


1 7+0 09 


1 8+0 03 


0.18+0.02 






1.50 


2.00 


0.25±0.06 


22+n 03 


n 1 ^+0 02 


1 5+0 03 

U . J- tJ—l—U . yJO 


0.15+0.02 






2.00 


2.50 


0.24±0.06 


0.23±0.03 


0.16+0.02 


0.15+0.03 


0.139+0.015 






2.50 


3.00 


± 


0.17±0.03 


0.15+0.02 


0.10+0.02 


0.127+0.012 






3.00 


3.50 


± 


0.12±0.02 


0.14+0.02 


0.14+0.02 


0.109+0.010 






3.50 


4.00 


± 


0.11+0.02 


0.094+0.011 


0.082+0.014 


0.086+0.008 






4.00 


5.00 


± 


0.11+0.02 


0.074+0.009 


0.052+0.009 


0.065+0.006 






5.00 


6.50 


± 


± 


0.036+0.005 


0.022+0.004 


0.027+0.003 






6.50 


8.00 


± 


± 


+ 


0.010+0.003 


0.011+0.002 


0.200 


0.250 


0.50 


1.00 


0.31±0.13 


0.26+0.06 


0.28+0.05 


0.17+0.05 


0.21+0.03 






1.00 


1.50 


0.12±0.05 


0.18+0.04 


0.20+0.04 


0.12+0.04 


0.16+0.02 






1.50 


2.00 


0.15±0.05 


0.14+0.03 


0.16+0.03 


0.11+0.03 


0.10+0.02 






2.00 


2.50 


0.10±0.04 


0.12+0.03 


0.14+0.03 


0.08+0.03 


0.09+0.02 






2.50 


3.00 


± 


0.11+0.02 


0.12+0.02 


0.06+0.02 


0.056+0.010 






3.00 


3.50 


± 


0.08+0.02 


0.064+0.011 


0.042+0.012 


0.039+0.006 






3.50 


4.00 


± 


0.050+0.011 


0.049+0.009 


0.032+0.010 


0.034+0.005 






4.00 


5.00 


± 


0.031+0.008 


0.045+0.008 


0.023+0.009 


0.030+0.005 






5.00 


6.50 


± 


± 


0.028+0.007 


0.006+0.004 


0.013+0.002 






6.50 


8.00 


± 


± 


± 


0.007+0.004 


0.007+0.001 
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TABLE XI: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, a'' / {dpdO,) , for 7r~-Cu interactions at 
3,5,8,12 GeV/c. Each row refers to a different {Pjjun < 
P < Pmaxj^'min < S < 6'max) bin, where p and 9 are 
the outgoing proton momentum and polar angle, re- 
spectively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 



e ■ 

mill 
(rad) 


Omax Pmin Pmax 
(rad) (GeV/c) (GeV/c) 


3 GeV/c 


tfa^/idpdn) 

(barn/(sr GeV/c)) 
5 GeV/c 8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.31±0.05 


0.35±0.05 


0.42±0.06 


0.33±0.07 






1.00 


1.50 


0.18±0.03 


0.29±0.03 


0.32±0.04 


0.31±0.05 






1.50 


2.00 


0.13±0.02 


0.21±0.03 


0.31±0.04 


0.26±0.05 






2.00 


2.50 


0.074±0.013 


0.14±0.02 


0.20±0.03 


0.09±0.03 






2.50 


3.00 


± 


0.075±0.013 


0.20±0.02 


0.17±0.03 






3.00 


3.50 


± 


0.033±0.008 


0.13±0.02 


0.22±0.03 






3.50 


4.00 


± 


0.021±0.006 


0.13±0.02 


0.17±0.03 






4.00 


5.00 


± 


0.016±0.005 


0.062±0.011 


0.15±0.02 






5.00 


6.50 


± 


± 


0.018±0.005 


0.071±0.010 






6.50 


8.00 


± 


± 


± 


0.034±0.007 


0.100 


0.150 


0.50 


1.00 


0.34±0.05 


0.45±0.07 


0.49±0.08 


0.45±0.09 






1.00 


1.50 


0.24±0.03 


0.35±0.04 


0.33±0.04 


0.41±0.06 






1.50 


2.00 


0.12±0.02 


0.25±0.03 


0.28±0.04 


0.21±0.04 






2.00 


2.50 


0.080±0.013 


0.17±0.02 


0.20±0.03 


0.12±0.03 






2.50 


3.00 


± 


0.068±0.012 


0.19±0.02 


0.12±0.03 






3.00 


3.50 


± 


0.050±0.009 


0.10±0.02 


0.16±0.03 






3.50 


4.00 


± 


0.016±0.007 


0.081±0.013 


0.12±0.02 






4.00 


5.00 


± 


0.011±0.004 


0.046±0.009 


0.09±0.02 






5.00 


6.50 


± 


± 


0.010±0.003 


0.026±0.006 






6.50 


8.00 


± 


± 


± 


0.013±0.004 


0.150 


0.200 


0.50 


1.00 


0.30±0.05 


0.38±0.06 


0.40±0.07 


0.36±0.08 






1.00 


1.50 


0.17±0.02 


0.20±0.03 


0.29±0.03 


0.25±0.04 






1.50 


2.00 


0.08±0.02 


0.16±0.02 


0.23±0.03 


0.19±0.04 






2.00 


2.50 


0.071±0.013 


0.13±0.02 


0.16±0.03 


0.17±0.04 






2.50 


3.00 


± 


0.063±0.011 


0.11±0.02 


0.12±0.02 






3.00 


3.50 


± 


0.033±0.008 


0.10±0.02 


0.12±0.02 






3.50 


4.00 


± 


0.009±0.004 


0.052±0.011 


0.08±0.02 






4.00 


5.00 


± 


0.003±0.002 


0.020±0.005 


0.041±0.010 






5.00 


6.50 


± 


± 


O.OOOiO.OOl 


0.009±0.003 






6.50 


8.00 


± 


± 


± 


0.003±0.002 


0.200 


0.250 


0.50 


1.00 


0.34±0.07 


0.32±0.07 


0.48±0.09 


0.38±0.09 






1.00 


1.50 


0.26±0.04 


0.33±0.05 


0.48±0.08 


0.45±0.09 






1.50 


2.00 


0.16±0.03 


0.32±0.05 


0.27±0.06 


0.29±0.07 






2.00 


2.50 


0.07±0.02 


0.19±0.03 


0.31±0.06 


0.31±0.08 






2.50 


3.00 


± 


0.06±0.02 


0.16±0.03 


0.19±0.04 






3.00 


3.50 


± 


0.024±0.008 


0.047±0.013 


0.09±0.02 






3.50 


4.00 


± 


0.006±0.003 


0.030±0.009 


0.039±0.012 






4.00 


5.00 


± 


0.003±0.002 


0.016±0.005 


0.017±0.007 






5.00 


6.50 


± 


± 


0.001±0.003 


0.004±0.003 






6.50 


8.00 


± 


± 


± 


0.002±0.002 



TABLE XII: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, (f(j^/{dpdQ), for TT^-Cu interactions at 
3,5,8,12 GeV/c. Each row refers to a different (Pjyiin < 
P < Pmax,6'nnn <9< 6'niax) bin, where p and 9 are 
the outgoing proton momentum and polar angle, re- 
spectively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 



e ■ 

mill 


0max 


Pmin 


pmax 




dPa'^/idpdO) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GcV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.25±0.12 




f) 43+n OS 


0.17+0.16 






1.00 


1.50 


0.28±0.10 




32+0 06 


0.48+0.21 






1.50 


2.00 


± 


25±0 05 


30+0 05 

yj — 1 — yj *yj^j 


0.30+0.41 






2.00 


2.50 


± 


n 05+0 05 


1 6+0 04 


0.19+0.32 






2.50 


3.00 


± 


± 


13+0 04 


0.12+0.16 






3.00 


3.50 


± 




OQ+0 04 

yj 'yjij 1 yj . 


0.07+0.07 






3.50 


4.00 


± 




03+0 04 


0.16+0.13 






4.00 


5.00 


± 


-I- 


± 


0.10+0.12 






5.00 


6.50 


± 


-1- 


-|- 


0.02+0.06 






6.50 


8.00 


± 


-I- 


+ 


0.01+0.07 


0.100 


0.150 


0.50 


1.00 


0.30±0.13 


n 47+ ry nq 


53+0 1 1 

yj - OO—I yj - J. J. 


0.69+0.36 






1.00 


1.50 


0.27±0.09 


\J .KjyJ -i—\J .\JiJ 




0.44+0.20 






1.50 


2.00 


0.07±0.08 


n 2fi+n 05 


0.24+0.05 


0.32+0.17 






2.00 


2.50 


± 


16+0 OQ 


0.14+0.04 


0.20+0.13 






2.50 


3.00 


± 


-I- 


0.13+0.03 


0.15+0.11 






3.00 


3.50 


± 


± 


0.12+0.03 


0.04+0.05 






3.50 


4.00 


± 


± 


0.04+0.02 


0.11+0.09 






4.00 


5.00 


± 


± 


0.04+0.01 


0.11+0.07 






5.00 


6.50 


± 




0.02+0.02 


0.04+0.04 






6.50 


8.00 


± 


-l- 


+ 


0.02+0.03 


0.150 


0.200 


0.50 


1.00 


0.27±0.12 


53+0 1 


0.46+0.10 


0.21+0.20 






1.00 


1.50 


0.27±0.08 


9S+0 04 


0.30+0.05 


0.34+0.18 






1.50 


2.00 


0.05±0.08 


n 1 Q+n 04 


0.20+0.04 


0.29+0.17 






2.00 


2.50 


± 


0.07±0.03 


0.15+0.04 


0.23+0.14 






2.50 


3.00 


± 


0.03±0.02 


0.12+0.03 


0.15+0.11 






3.00 


3.50 


± 


0.01±0.02 


0.10+0.03 


0.08+0.09 






3.50 


4.00 


± 


0.02+0.02 


0.05+0.02 


0.11+0.11 






4.00 


5.00 


± 


0.01+0.02 


0.02+0.01 


0.02+0.04 






5.00 


6.50 


± 


± 


0.003+0.003 0.001+0.003 






6.50 


8.00 


± 


+ 


+ 


0.006+0.025 


0.200 


0.250 


0.50 


1.00 


0.32±0.17 


0.33+0.10 


0.39+0.10 


0.05+0.12 






1.00 


1.50 


0.17±0.09 


0.30+0.06 


0.44+0.09 


0.09+0.17 






1.50 


2.00 


0.08±0.05 


0.12+0.04 


0.16+0.05 


0.04+0.14 






2.00 


2.50 


0.07±0.07 


0.10+0.04 


0.14+0.05 


0.26+0.31 






2.50 


3.00 


± 


0.06+0.02 


0.09+0.03 


0.09+0.22 






3.00 


3.50 


± 


0.04+0.02 


0.04+0.03 


0.01+0.11 






3.50 


4.00 


± 


0.03+0.01 


0.02+0.02 


+ 






4.00 


5.00 


± 


0.02+0.01 


0.024+0.035 


+ 






5.00 


6.50 


± 


± 


0.009+0.012 


+ 






6.50 


8.00 


± 


± 


± 


+ 



TABLE XIII: HARP results for the double-differential 
p production cross-section in the laboratory system, 
(f a'' / (dpdQ,) , for p-Cu interactions at 3,5,8,12 GeV/c. 
Each row refers to a different (Pjnin < P < 
Pmax,^niin ^ < ^max) bin, where p and 9 arc the 
outgoing proton momentum and polar angle, respec- 
tively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 



e ■ 

mill 


Omax 


Pmin 


pmax 




d^cr^/idpcl^}) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GcV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.48±0.19 


n 4.5+n nq 


f) 53+0 ns 


0.51±0.09 






1.00 


1.50 


0.36±0.11 


n 49+0 06 


59+0 05 


0.54+0.07 






1.50 


2.00 


0.40±0.11 


54±0 07 


54+0 05 


0.46+0.06 






2.00 


2.50 


0.46±0.15 


57+0 Ofi 


50+0 04 


0.44+0.05 






2.50 


3.00 


± 


76±0 08 


55+0 05 


0.57±0.05 






3.00 


3.50 


± 


n 87+n 07 


64+0 04 


0.49+0.05 






3.50 


4.00 


± 


77+0 07 


66+0 05 


0.52+0.05 






4.00 


5.00 


± 


f) sfi+n 1 2 


n sn+n 04 


0.62±0.04 






5.00 


6.50 


± 


-1- 


n c;4_(_n no 


0.51+0.03 






6.50 


8.00 


± 


-I- 


± 


0.41+0.02 


0.100 


0.150 


0.50 


1.00 


0.40±0.17 


n 54.+n 1 n 


67+0 OQ 


0.69±0.12 






1.00 


1.50 


0.59±0.15 


n 4Q-|-n (17 


48+0 05 


0.48±0.07 






1.50 


2.00 


0.62±0.15 


n 45+n 06 


0.47+0.05 


0.43+0.06 






2.00 


2.50 


0.52±0.12 


50+n 07 


0.41+0.04 


0.44+0.06 






2.50 


3.00 


± 


52+0 06 


0.43±0.05 


0.39±0.05 






3.00 


3.50 


± 


50+0 06 


0.49+0.04 


0.35+0.04 






3.50 


4.00 


± 


45+0 05 


0.42+0.03 


0.40+0.04 






4.00 


5.00 


± 


50+0 06 


0.34±0.03 


0.35±0.03 






5.00 


6.50 


± 




0.19±0.02 


0.20±0.02 






6.50 


8.00 


± 


+ 


+ 


0.104+0.013 


0.150 


0.200 


0.50 


1.00 


0.37±0.16 


54+0 1 1 


0.61+0.10 


0.57+0.11 






1.00 


1.50 


0.51±0.13 


4n+0 06 

\J .^yj^^ZyJ .\J\J 


0.40±0.05 


0.35±0.05 






1.50 


2.00 


0.45±0.13 


4'^+0 05 


0.32+0.04 


0.29+0.05 






2.00 


2.50 


0.36±0.11 


0.42±0.06 


0.38+0.04 


0.29+0.05 






2.50 


3.00 


± 


0.23±0.04 


0.28±0.03 


0.17±0.03 






3.00 


3.50 


± 


0.16±0.03 


0.21±0.03 


0.22±0.03 






3.50 


4.00 


± 


0.17±0.03 


0.16+0.02 


0.16±0.02 






4.00 


5.00 


± 


0.16+0.02 


0.14+0.02 


0.12+0.02 






5.00 


6.50 


± 


± 


0.051±0.008 0.034±0.007 






6.50 


8.00 


± 


+ 


+ 


0.016+0.004 


0.200 


0.250 


0.50 


1.00 


0.65±0.29 


0.44+0.11 


0.59+0.10 


0.62+0.13 






1.00 


1.50 


0.23±0.12 


0.43±0.08 


0.40±0.07 


0.36±0.08 






1.50 


2.00 


0.21±0.09 


0.29+0.06 


0.32+0.06 


0.20+0.06 






2.00 


2.50 


0.27±0.11 


0.18+0.04 


0.26+0.05 


0.18+0.05 






2.50 


3.00 


± 


0.19+0.04 


0.18+0.03 


0.12+0.04 






3.00 


3.50 


± 


0.08±0.02 


0.13±0.02 


0.09±0.02 






3.50 


4.00 


± 


0.07+0.02 


0.067+0.014 


0.05+0.02 






4.00 


5.00 


± 


0.06+0.02 


0.063+0.013 


0.043+0.014 






5.00 


6.50 


± 


± 


0.024±0.006 0.027±0.010 






6.50 


8.00 


± 


± 


± 


0.019±0.008 
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TABLE XIV: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, dF / [dpdO) , for 7r~-Sn interactions at 
3,5,8,12 GeV/c. Each row refers to a different 
(Pmin < -P < Pmax,^'rnin < 9 < ^'max) bin, where p 
and 9 are the outgoing proton momentum and polar 
angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance 
matrix are given. 



9 ■ 
mill 

(rad) 


0max Prnin Pmax 
(rad) (GeV/c) (GeV/c) 


3 GeV/c 


(fa^/idpdO) 
(barn/(sr GeV/c)) 
5 GeV/c 8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.52±0.09 


0.67±0.10 


0.79±0.11 


0.75±0.11 






1.00 


1.50 


0.26±0.05 


0.45±0.05 


0.51±0.07 


0.55±0.06 






1.50 


2.00 


0.15±0.03 


0.30±0.04 


0.47±0.06 


0.41±0.08 






2.00 


2.50 


0.05±0.02 


0.19±0.03 


0.24±0.05 


0.11±0.34 






2.50 


3.00 


± 


0.09±0.02 


0.26±0.03 


0.27±0.07 






3.00 


3.50 


± 


0.043±0.013 


0.12±0.03 


0.33±0.07 






3.50 


4.00 


± 


0.029±0.009 


0.14±0.02 


0.25±0.04 






4.00 


5.00 


± 


0.014±0.007 


0.07±0.02 


0.22±0.03 






5.00 


6.50 


± 


± 


0.019±0.007 


0.121±0.015 






6.50 


8.00 


± 


± 


± 


0.046±0.010 


0.100 


0.150 


0.50 


1.00 


0.53±0.09 


0.62±0.10 


0.82±0.13 


0.81±0.13 






1.00 


1.50 


0.33±0.06 


0.55±0.07 


0.61±0.08 


0.72±0.08 






1.50 


2.00 


0.16±0.03 


0.32±0.04 


0.38±0.06 


0.42±0.07 






2.00 


2.50 


0.06±0.02 


0.21±0.03 


0.32±0.04 


0.22±0.07 






2.50 


3.00 


± 


0.10±0.02 


0.21±0.03 


0.30±0.04 






3.00 


3.50 


± 


0.050±0.012 


0.15±0.03 


0.22±0.04 






3.50 


4.00 


± 


0.019±0.007 


0.12±0.02 


0.17±0.02 






4.00 


5.00 


± 


0.017±0.007 


0.041±0.011 


0.14±0.02 






5.00 


6.50 


± 


± 


0.009±0.003 


0.052±0.010 






6.50 


8.00 


± 


± 


± 


0.011±0.003 


0.150 


0.200 


0.50 


1.00 


0.51±0.10 


0.82±0.13 


0.66±0.12 


0.58±0.11 






1.00 


1.50 


0.19±0.04 


0.37±0.05 


0.44±0.06 


0.58±0.07 






1.50 


2.00 


0.09±0.03 


0.28±0.04 


0.37±0.05 


0.24±0.05 






2.00 


2.50 


0.04±0.02 


0.15±0.03 


0.21±0.05 


0.22±0.10 






2.50 


3.00 


± 


0.07±0.02 


0.17±0.03 


0.27±0.07 






3.00 


3.50 


± 


0.024±0.008 


0.13±0.03 


0.15±0.03 






3.50 


4.00 


± 


0.007±0.004 


0.08±0.02 


0.09±0.02 






4.00 


5.00 


± 


0.003±0.003 


0.020±0.006 


0.071±0.014 






5.00 


6.50 


± 


± 


O.OOliO.OOl 


0.011±0.004 






6.50 


8.00 


± 


± 


± 


0.002±0.001 


0.200 


0.250 


0.50 


1.00 


0.57±0.13 


0.69±0.13 


0.83±0.16 


0.73±0.14 






1.00 


1.50 


0.29±0.07 


0.60±0.09 


0.86±0.14 


0.70±0.12 






1.50 


2.00 


0.14±0.04 


0.35±0.06 


0.48±0.10 


0.53±0.11 






2.00 


2.50 


0.07±0.03 


0.18±0.04 


0.43±0.08 


0.37±0.10 






2.50 


3.00 


± 


0.11±0.03 


0.17±0.05 


0.27±0.06 






3.00 


3.50 


± 


0.025±0.012 


0.04±0.02 


0.12±0.04 






3.50 


4.00 


± 


0.005±0.003 


0.016±0.009 


0.06±0.02 






4.00 


5.00 


± 


0.003±0.002 


0.007±0.004 


0.04±0.02 






5.00 


6.50 


± 


± 


0.001±0.006 


0.008±0.004 






6.50 


8.00 


± 


± 


± 


0.002±0.002 



TABLE XV: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, (P a'' / {dpdO,) , for Tr^-Sn interactions at 
3,5,8,12 GeV/c. Each row refers to a different 
(Pmin < -P < Pmax,^'rnin < 8 < dm&x) bin, where p 
and 9 are the outgoing proton momentum and polar 
angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance 
matrix are given. 



e ■ 

mill 


Omax 


Pmin 


pmax 








(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GcV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.49±0.16 


n 7n-|-n 1 s 


f) fifi+n 1 2 


0.75+0.27 






1.00 


1.50 


0.23±0.10 


n "14+0 07 


69+0 OQ 


0.54+0.17 






1.50 


2.00 


± 


37±0 12 


50+0 08 


0.45+0.60 






2.00 


2.50 


± 


n 01 +0 03 


36+0 07 


0.21+0.32 






2.50 


3.00 


± 


± 


15±0 06 

\J • ^^\J m\J\J 


0.26+0.25 






3.00 


3.50 


± 


± 


05+0 1 1 


0.21+0.11 






3.50 


4.00 


± 


± 


01+0 33 


0.21+0.10 






4.00 


5.00 


± 


± 


-I- 


0.12+0.08 






5.00 


6.50 


± 




-|- 


+ 






6.50 


8.00 


± 


-I- 


-I- 


+ 


0.100 


0.150 


0.50 


1.00 


0.55±0.17 


n H4.+n 1 5 


1 03+0 1Q 


0.88+0.30 






1.00 


1.50 


0.43±0.12 


n CI -(.(-) no 




0.76+0.21 






1.50 


2.00 


0.14±0.12 


n 41 +n 08 


0.48+0.08 


0.36+0.14 






2.00 


2.50 


± 


10+0 06 


0.22+0.05 


0.33+0.14 






2.50 


3.00 


± 


± 


0.20±0.04 


0.19+0.10 






3.00 


3.50 


± 


± 


0.21+0.06 


0.08+0.06 






3.50 


4.00 


± 


± 


0.08+0.04 


0.15+0.09 






4.00 


5.00 


± 


± 


0.08+0.04 


0.14+0.07 






5.00 


6.50 


± 


-|- 


0.01+0.02 


0.05+0.06 






6.50 


8.00 


± 


-|- 


+ 


0.02+0.04 


0.150 


0.200 


0.50 


1.00 


0.64±0.19 


Qfi+0 16 


0.67+0.14 


0.74+0.30 






1.00 


1.50 


0.29±0.08 


n 4.9-l-n Of, 


0.48+0.08 


0.63+0.19 






1.50 


2.00 


± 


n 30+n 05 

\j .fjyj -i—U .yjo 


0.38+0.07 


0.43+0.16 






2.00 


2.50 


± 


0.10±0.05 


0.23+0.05 


0.25+0.12 






2.50 


3.00 


± 
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TABLE XVI: HARP results for the double-differential 
p production cross-section in the laboratory system, 
<fa^/{dpdQ), for p-C interactions at 3,5,8,12 GeV/c. 
Each row refers to a different (Pj^i^ < P < 
Pmaxj^'inin < ^ < ^max) bin, whore p and 6 arc the 
outgoing proton momentum and polar angle, respec- 
tively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 
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TABLE XVII: HARP results for the double- 
differential p production cross-section in the labora- 
tory system, cf a'^ / (dpdQ) , for 7r~-Ta interactions at 
3,5,8,12 GeV/c. Each row refers to a different {Pjjun < 
P < Pmaxj^'min < S < 6'max) bin, where p and 9 are 
the outgoing proton momentum and polar angle, re- 
spectively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 
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TABLE XVIII: HARP results for the double- 
differential p production cross-section in the labora- 
tory system, cP / [dpdQ) , for Tr^-Ta interactions at 
3,5,8,12 GeV/c. Each row refers to a different (Pjyiin < 
P < Pmaxj^'min < S < 6'max) bin, where p and 9 are 
the outgoing proton momentum and polar angle, re- 
spectively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 
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TABLE XIX: HARP results for the double-differential 
p production cross-section in the laboratory system, 
<fa^/{dpdQ), for p-Ta interactions at 3,5,8,12 GeV/c. 
Each row refers to a different (Pj^i^ < P < 
Pmaxj^'inin < ^ < ^max) bin, whore p and 6 arc the 
outgoing proton momentum and polar angle, respec- 
tively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 
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TABLE XX: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, (f (j^ / {dpdO) , for 7r~-Pb interactions at 
3,5,8,12 GeV/c. Each row refers to a different (Pjyiin < 
P < Pmaxj^'min < S < dinux) bin, where p and 9 are 
the outgoing proton momentum and polar angle, re- 
spectively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 
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TABLE XXI: HARP results for the double-differential 
p production cross-section in the laboratory sys- 
tem, (P a'' / (dpdQ) , for Tr^-Pb interactions at 
3,5,8,12 GeV/c. Each row refers to a different {Pjjun < 
P < Pmaxj^'min < S < 6'max) bin, where p and 9 are 
the outgoing proton momentum and polar angle, re- 
spectively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 
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TABLE XXII: HARP results for the double-differential 
p production cross-section in the laboratory system, 
<fa^/{dpdQ), for p-Pb interactions at 3,5,8,12 GeV/c. 
Each row refers to a different (Pj^i^ < P < 
Pmaxj^'inin < ^ < ^max) bin, whore p and 6 arc the 
outgoing proton momentum and polar angle, respec- 
tively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are 
given. 
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mill 


0max 


Pmin 


pmax 








(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GeV/c)) 
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1 02+0 09 


1.07+0.15 






3.50 


4.00 


± 


1 05+0 1 5 


1 12+0 08 


1.09+0.12 






4.00 


5.00 


± 


f) sfi+n 1 1 

W.U"J_I_W. -L J- 


1 04+0 07 


0.97+0.09 






5.00 


6.50 


± 


-|- 


7Q+n 05 


0.83+0.07 






6.50 


8.00 


± 


-I- 


+ 


0.63+0.06 


0.100 


0.150 


0.50 


1.00 


0.59±0.20 


1 40-1-0 25 


1 QI +0 25 


2.39+0.40 






1.00 


1.50 


0.87±0.20 




1.17+0.13 


1.56+0.22 






1.50 


2.00 


0.45±0.15 


84+n 1 3 


92+0 1 1 


1.26+0.19 






2.00 


2.50 


0.39±0.10 


77+0 1 5 


70+0 08 


0.59+0.11 






2.50 


3.00 


± 


f) qy+n 1 3 

\J 'Cf 1 _1_ \J • -LlJ 


7«+0 08 


0.66+0.11 






3.00 


3.50 


± 


n 87+n 12 


67+0 07 


0.61+0.10 






3.50 


4.00 


± 


fifi+0 1 3 


60+0 06 


0.70+0.10 






4.00 


5.00 


± 


n t^q+n 1 n 


n 55+n 05 


0.54+0.07 






5.00 


6.50 


± 


-|- 


9Q+0 03 


0.35+0.04 






6.50 


8.00 


± 


-|- 


±- 


0.17+0.03 


0.150 


0.200 


0.50 


1.00 


0.93±0.30 


1 52+0 28 


1 55+0 23 


1.83+0.34 






1.00 


1.50 


0.52±0.14 




q4+o 1 

VJ. J'±^tVJ. J.\j 


1.04+0.17 






1.50 


2.00 


0.21±0.09 


n 63+0 1 1 


65+0 OS 


0.79+0.15 






2.00 


2.50 


0.37±0.13 


0.58±0.11 


0.50+0.07 


0.76+0.14 






2.50 


3.00 


± 


0.42±0.08 


0.43+0.06 


0.48+0.09 






3.00 


3.50 


± 


0.31±0.06 


0.30+0.04 


0.44+0.08 






3.50 


4.00 


± 


0.33+0.06 


0.22+0.03 


0.40+0.07 






4.00 


5.00 


± 


0.22+0.04 


0.23+0.03 


0.18+0.04 






5.00 


6.50 


± 


± 


0.058+0.013 


0.13+0.03 






6.50 


8.00 


± 


+ 


± 


0.06+0.02 


0.200 


0.250 


0.50 


1.00 


0.50±0.30 


1.49+0.33 


1.56+0.25 


2.05+0.43 






1.00 


1.50 


0.30±0.14 


0.84+0.18 


0.89+0.15 


0.89+0.22 






1.50 


2.00 


0.21±0.08 


0.30+0.10 


0.49+0.10 


0.63+0.18 






2.00 


2.50 


0.13±0.06 


0.30+0.09 


0.37+0.08 


0.41+0.14 






2.50 


3.00 


± 


0.35+0.09 


0.35+0.06 


0.31+0.10 






3.00 


3.50 


± 


0.20+0.06 


0.19+0.04 


0.17+0.06 






3.50 


4.00 


± 


0.09+0.03 


0.12+0.03 


0.11+0.04 






4.00 


5.00 


± 


0.10+0.03 


0.11+0.03 


0.13+0.05 






5.00 


6.50 


± 


± 


0.038+0.015 


0.09+0.04 






6.50 


8.00 


± 


± 


± 


0.05+0.02 



